A theory for magnetic-field effects of nonmagnetic organic semiconducting materials 
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A universal mechanism for strong magnetic-field effects of nonmagnetic organic semiconductors 
is presented. A wealc magnetic field (less than hundreds mT) can substantially change the charge 
carrier hopping coefficient between two neighboring organic molecules when the two hopping states 
are not too symmetric. Under the illumination of lights or under a high electric field, the change 
of hopping coefficients leads also to the change of polaron density so that photocurrent, photolumi- 
nescence, electroluminescence, magnetoresistance and electrical-injection current become sensitive 
to a weak magnetic field. The present theory can not only explain all observed features, but also 
provide a solid theoretical basis for the widely used empirical fitting formulas. 



One of the long-term pj-lJ] unsolved fundamental issues 
in organic physics is the mechanism behind the strong 
responses of electrical and optical properties of nonmag- 
netic organic semiconductors to a weak magnetic field, 
known as organic magnetic-field effect (OMFE). The re- 
cent revival interest in OMFE of organic semiconductors 
is largely due to its importance in fundamental science 
and technology applications Firstly, there is a be- 
lief that the OMFE can be used as a powerful tool to 
probe microscopic processes of organic materials. Sec- 
ondly, the OMFE can be used to develop new multi- 
functional organic devices f^. Experiments showed that 
OMFE has following surprising yet universal features. 

1) The OMFE appears in vast different organic semi- 
conductors without any magnetic elements at room tem- 
perature although the possible energy level shifts due to 
the presence of a magnetic field are orders magnitude 
smaller than the thermal energy and other energy scales. 

2) The electroluminescence, photocurrent, photolumines- 
cence, and electrical-injection current are very sensitive 
to weak magnetic field with both positive and negative 
OMFE though positive OMFE (or negative magnetore- 
sistance (MR) in the convention terminology) at very 
weak field is typically observed. 3) The OMFE can of- 
ten be fitted by two empirical formulas: [B/{B -I- Bq)Y 
and /{B^ + Bq)[Q], where B is the applied magnetic 
field. In the theoretical side, it is knownf?, ^ that famil- 
iar MR mechanisms such as Lorentz force, conventional 
hopping MR, electron-electron interaction and weak lo- 
calization are highly unlikely to be the cause behind the 
OMFE. The current belief in the community is that the 
OMFE is intimately tied to spin physics involving spin 
configuration, spin correlation, and spin flip 4]. How- 
ever, there is no convincing arguments why an extremely 
small Zeeman energy can beat other much larger energy 
scales in controlling electron spin dynamics to generate 
this OMFE. Both extensive experimental and theoreti- 
cal studies so far are suggesting that a novel explanation 
is needed. This new MR mechanism should explain not 



only all OMFE features, but also why the similar effects 
do not often appear in the usual inorganic semiconduc- 
tors. In this report, we present such a theory that does 
not explicitly rely on the electron spin degrees of free- 
dom. It is showed that the OMFE originates from the 
substantial change of electron hopping coefficient in a 
magnetic field because of narrow bandwidth of organic 
semiconductors and asymmetry in organic molecules. 

Organic semiconductors have a few distinct properties 
that their inorganic counterparts do not have. Firstly, 
unlike an atom that is sphere-like, an organic molecule 
is highly irregular. As a result, organic molecular wave- 
function has no obvious symmetry. Secondly, organic 
molecules in organic semiconductors are bonded by the 
Van der Waals force so that their bands are very narrow 
in comparison with an order of lOeV bandwidth for their 
inorganic counterparts [ll| . Thirdly, the intramolecular 
excitons have strong binding energies of order of eV[4|. 
On the other hand, the electron and hole become po- 
laron pair when they are located on different molecules 
because of weak intermolecular exciton binding energy Q . 
The electrical properties of an organic semiconductor are 
mainly determined by the motion of polarons since the 
motion of excitons does not contribute to the electric 
current. The singlet excitons are responsible to the lu- 
minescence. A weak field should not change much of 
energy levels of various excitation states so that their 
populations at thermal equilibrium are not sensitive to a 
magnetic field because they depend only on the energy 
level distribution and the temperature. Any significant 
change in magnetoresistance near the quasiequilibrium 
state must come from the mobility change. The question 
is whether a weak field of lOOmT can change the mobil- 
ity of polarons in organic semiconductors. In the usual 
inorganic crystals with s-like wavefunction, the answer is 
no. However, we argue below that this can indeed hap- 
pen in organic conjugated materials with highly irregular 
molecular wavefunction. 

In order to understand why a weak magnetic field can 
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FIG. 1: Schematical draw of two organic molecules (Alqa) sep- 
arated by a distance d and aligned along x-direction. Molecule 
1 is centered at the origin, and molecule 2 is centered at 
(d, 0,0). The field is assumed to be along z-direction. tpi and 
^2 are two localized states with localization lengths and ^2 
on molecules 1 and 2, respectively. 



change charge carrier (electron and hole or polaron) mo- 
bility in an organic conjugated material, we consider a 
system with two molecules separated by a distance d as 
schematically shown in Fig. 1. One-electron Hamiltonian 
in a magnetic field can in general be described by 



2m c 



V1+V2 



(1) 



where Vi and V2 are the potential created by molecules 
1 and 2, respectively. A is the vector potential due 
to magnetic field B. For the simplicity and clarity, we 
shall assume that the two molecules are aligned along 
x-direction, the field is along the z-direction (pointing 
out of the paper). The important quantity for electron 
transport is the tunneling matrix element between two 
molecules. When an electron tunnels from an initially 
occupied state, say ipi of molecule 1, to empty state 
■02 of molecule 2 with tunneling matrix t , it will con- 
tribute to the hopping probability P (per unit time), 
proportional to exp(— Aei2/(i4rT)), where Aei2 de- 
scribes the relative energy level with respect to the Fermi 
level [12l|. The hopping conduction can be regarded as 
an electron diffusion process with a diffusion constant 
D = Pd^, where d is the average distance between two 
neighboring molecules. According to the Einstein rela- 
tion, the electron mobility /i is given hy fi — eD/{KT) 
which is related to the conductivity in the conventional 
wav[l^. Therefore, we can concentrate on how the tun- 
neling matrix element depends on the magnetic field in 
order to study the magnetoresistance of the system. 

In the tight-binding approximation ^3|, one of the au- 
thors in an early publication has generalized the 
Bardeen's transfer matrix formalism to high dimension 
and in the presence of a magnetic field. In 3D, it is 

(2) 

where 4>o = ch/e is the flux quanta. For small A when 
the magnetic length Ib = yVo/^ is bigger than d, mag- 
netic confinement that is responsible for the exponential 



increase of resistance in the usual hopping conduction 
can be neglected and ■01 and ip2 do not depend on B to 
the zero order approximation. Then the magnitude of 
the field-independent part of t is order of 
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J '^*iMx=idydz 



while that of the field dependent part is 

^TT I y^i^^L^idydz. 
m J ^ 

= -I- ^2~^, and ^1 and ^2 are the localization 
lengths of "01 and ■02 , respectively. Both terms depend on 
the nature of the wavefunctions. This explains why the 
OMFE value varys from sample to sample, and from ma- 
terial to material ||4i]. The second term vanishes for s-like 
wavefunction. This is why the similar phenomena do not 
show up in usual inorganic semiconductors. Due to the 
irregular structures of organic molecules, one will expect 
an appreciable value for the second term, resulting in a 
sizable change of the hopping probability. 

Due to the Van der Waals bonding, the OMFE is 
measurable only under an optical injection of carriers or 
an electric carrier injection by an electric field above a 
threshold. When an organic semiconductor is under the 
illumination of a light or under a high electric field, the 
field dependent t results in a field dependence of polaron 
density. Take optical injection of carriers as an example, 
under the illumination of a light, an electron in a high- 
est occupied molecular orbit (HOMO) absorbs a photon 
and jumps to a higher empty molecular orbit of the same 
molecule. As schematically illustrated in Fig. 2, the ex- 
cited electron can either dump its excessive kinetic energy 
to its environment and forms an exciton with the hole left 
behind or jumps to neighboring molecules and becomes 
polarons. Depending on the relative probabilities of ex- 
cited electrons (holes) staying in the same molecules and 
jumping to different molecules, the polaron density shall 
vary with the illumination intensity. Let us denote the 
probability (per unit time) of a pair of electron and hole 
on the same molecule forming an exciton by Pq ~ ^/t, 
where r is the typical time for a pair of electron and 
hole to form an exciton. Pq is not sensitive to a weak 
field since the field cannot change much molecule orbits 
that determine Pq. Then the polaron generation rate 
per unit volume is JP/(Po -I- P) where J is the pho- 
ton absorption rate per unit volume and P cx |ip is the 
intermolecular hopping probability. Without the illumi- 
nation of a light, polaron density shall reach its equilib- 
rium density tiq at a rate of 7(71 — tiq), where 7 is po- 
laron decay rate. At balance, JPj (Pq -I- P) = 7(71 — no), 
thus the photon-generated polaron density n should be 
no -I- JPj [7(Po -I- P)] . Clearly, P-dependence of P results 
in a P— dependence of polaron density. 

Since the OMFE can often be fitted by two empirical 
functions B'^/{B'^ + B^) and [B/{B + Po)]^L6j|, a correct 
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FIG. 2: Schematic illustration of polaron and exciton forma- 
tion after a pair of electron and hole is created by a photon 
absorption. The excited electron-hole pair has probability 
P jumping to the neighboring molecules to form positive- 
charged and negative-charged polarons, and probability Po 
to form an exciton. 



theory should be able to provide a theortical basis for the 
fact. According to Eq. ([2]), t takes a form of Bq + iaB 
with Bq and a real and field-independent parameters if 
ipi and V'2 are real functions. This is the case when the 
molecule orbits involved in hopping are localized or not 
degeneratedfis'l. In this case, P oc |ip = a{B'^+Bl) and 
the polaron density shall depend on the magnetic field as 
^{Po+p) J + "0 = + aB^/{B^ + Bl), where no, nf,, a 
and Bq are B-independent parameters that depend on the 
molecule orbits involved. Thus, B'^ / {B^ + Bq) is a natu- 
ral OMFE function for t = Bo + iaB. The second type of 
empirical function appears naturally for t — i^Bq + aB). 
According to Eq. ([2]), this can happen when the spa- 
tial derivatives of tpi or ^2 are the functions multiplied 
by pure imaginary numbers. Of course, this must cor- 
respond to degenerated states. In this case, the leading 
term in the polaron density takes a form of [B / {B + Bq)Y 
in a similar argument when P ^ Pq. In reality, electron 
(polaron) hopping between two organic molecules should 
involve many molecule orbits, especially in photophysical 
processes and in a high electric field. One then needs to 
add contributions from all hopping events. Thus, it is 
likely that both B^/{B^ + B^) and [B/{B + Bq)]^ pro- 
cesses are presented, and OMFE should then be fitted by 
the linear combinations of these two functions, consistent 
with experimental findings. 

The novel mechanism is very robust. At the room 
temperature, the transport of charge carriers will involve 
many different molecule orbits. Each hopping event will 
subject to the influence of this mechanism as long as 
magnetic confinement is negligible {Ib > d) and molecu- 
lar structure is not sphere-like. Of course, thermal aver- 
age over all hopping events is needed. Molecule-molecule 
orientation in organic semiconductors should be quasi- 
random due to the nature of organic molecules. A mag- 
netic field can be along any direction with respect to the 



molecule-molecule bond instead of perpendicular direc- 
tion as assumed in the above discussion. This explains 
why OMFE is not sensitive to the field direction in de- 
vices. According to Eq. ([2]), different angle between 
the field and molecule-molecule bond leads to different 
hopping coefficient. It should also be emphasized that 
the mechanism present here does not depend on electron 
spins, and it does not require large energy splits of differ- 
ent spin configurations. It is applicable to both bipolar 
and hole-only (or electron-only) devices. Differ from the 
previous theories that try to relate the OMFE to the 
changes of electron levels, the present theory attributes 
the OMFE to the change of electron hopping coefficient 
in a field. Thus, it does not have all the troubles as 
those spin-dynamics related theories involve concepts of 
excitons and bipolarons 

In conclusion, we present a novel mechanism for the 
OMFE for nonmagnetic organic semiconductors. The 
mechanism is very general and robust for organic semi- 
conductors, but is normally not important for usual cova- 
lently bonded inorganic semiconductors. The mechanism 
can not only explain all experimentally observed OMFE, 
but also provides a solid theoretical basis for the empir- 
ical OMFE formulas. New experiments are needed to 
firmly establish this mechanism as the genuine cause of 
the OMFE. 
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